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ABSTRACT 

We analyze two regions of the quiet Sun (35.6 x 35.6 Mm 2 ) observed at high spatial resolution 
(<100 km) in polarized light by the IMaX spectropolarimeter onboard the Sunrise balloon. We identify 
497 small-scale (~400 km) magnetic loops, appearing at an effective rate of 0.25 loop h _1 arcsec -2 ; 
further, we argue that this number and rate are underestimated by ^30%. However, we find that 
these small dipoles do not appear uniformly on the solar surface: their spatial distribution is rather 
filamentary and clumpy, creating dead calm areas, characterized by a very low magnetic signal and 
a lack of organized loop-like structures at the detection level of our instruments, that cannot be 
explained as just statistical fluctuations of a Poisson spatial process. We argue that this is an intrinsic 
characteristic of the mechanism that generates the magnetic fields in the very quiet Sun. The spatio- 
temporal coherences and the clumpy structure of the phenomenon suggest a recurrent, intermittent 
mechanism for the generation of magnetic fields in the quietest areas of the Sun. 
Subject headings: Sun: surface magnetism — Sun: dynamo — Polarization 



1. INTRODUCTION 

During the last few years, our understanding of 
the structure, organization and evolution of mag- 
netic fields in the very quiet Sun (the regions out- 
side active regions and the network) has become in- 
creasingly clear. Magnetic fields in the quietest ar- 
eas of the Sun are relatively weak and organized at 
small spatial scales, which yields weak polarization 
signals that are difficult to observe. Until very re- 
cently, the general picture of the structure of its mag- 
netis m was rather r ough: a "turbulent" disorganized 
field (IStenflol IT982T: ISolankil 119931: IManso Sainz et al.l 
120041: iTruiillo Bueno et all 120041 ). It is now clear that 
even in very quiet areas, magnetic fields may organize 
as coherent loops at granular and subgra nular scales 
(<1000 km; IMartmez Gonzalez et all 120071). that these 
small loops aredynamic ([Martinez Gonzalez et alj|2007t 
Centeno et al.l I2007t IMa r tmez Gonzalez fc Bellot Rubiol 
20091: IGomorv et alJl2010T) . that thev pervade the quiet 
solar surface a nd may even connect with upper atmo- 
spheric layers (IMartmez Gonzalez fc Bellot Rubiol 120091 : 
IMartmez Gonzalez et al.ll2010D .~ 

Yet, this picture is still incomplete. For example, we 
lack a complete mapping of the full magnetic field vector 
on extended fields of view, because the linear polariza- 
tion signals are intrinsically weak (they are second order 
on the transverse magnetic field component), and high 
spatial resolution maps on lin ear polarization are rather 
patchy (Danil ovic et al.ll2010D . which has led to incom- 
plete (and sometimes, physically pr oblematic) character- 
izations of the topology of the fi eld (|Ishikawa et al.ll2008l : 
llshikawa fc Tsimetdl200lll20Tol) . 

Here we look for and trace small magnetic loops on ex- 
tended regions of the quiet Sun observed with the high- 
est spatial resolution. Loop-like structures are a natural 



configuration of the magnetic field due to its solenoidal 
character. While they can be traced as single, individual, 
coherent entities, they characterize the magnetic field at 
large, and their statistics and evolution may shed some 
light on the origin of the very quiet Sun magnetism, in 
parti cular, on the o peration or not of local dynamo ac- 
tion (jCatta nco 1999). On the other hand, the organiza- 
tion of the field at small scales affect the organization of 
the magnetic fi eld at larger scales and in higher atmo - 
spherics layers (iSchriiver fc Title] 120031: iCranmerl 120091) , 
and dynamics (|Cranmer fc van Ballegooiiej2^10l ). 

We find evidence for the small scale loops appearing 
rather irregularly, as in bursts and clumps. Moreover, 
wide regions of the very quiet Sun show very low mag- 
netic activity and no apparent sign of organized loops at 
the detection level of the instruments. These extremely 
quiet (dead calm) regions are an intrinsic characteristic 
of the statistical distribution of these events. 

2. IMAX DATA 

This paper is devoted to the analysis of disk cen- 
ter quiet Sun observations obtained with the IMaX in- 
strument (|Martmez Pillet et al.ll2011l) onboard the SUN- 
RISE balloon bor ne observatory ()Solanki et all 120101 : 
IBarthol et "aIll2010D . IMaX is a Fabry-Perot interferome- 
ter with polarimetric capabilities at the Fei line at 5250.2 
A. We analyze two different data sets. Both have the 
same properties except that they trace different regions 
of the quiet Sun and were observed at different times, 
being times series of 22 and 31 min duration. They con- 
sist of five filtergrams taken at ±40, ±80 and +227 mA 
from the Fei 5250.2 A line center. The field of view 
is 46.8" x 46.8" (35. 6x35.6 Mm 2 ), 20 times larger than 
the on e observed by IMartmez Gonzalez fc Bellot Rubiol 
(2009). The spatial resolution is of about 0".15- 
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Fig. 1. — Four examples of small loops found in the IMaX data. 
The grayscale background images represent the local granulation. 
The red and blue contours represent isocontours of magnetic flux 
density delimiting the positive and negative footpoints of the loop, 
respectively. The yellow contours represent the linear polarization. 
The time goes from top to bottom and the time intervals between 
successive images are not regular. The values of the polarization 
contours are different but are above the noise level and are chosen 
in each case to clearly see the loop. 

0" .18". The time cadence is 32 s (note th at it is 28 s 
in iMartinez Go nzalez & Bcllot Rubio 2009), allowing a 
noise level of a — 10 -3 and 7 x 10~ 4 I c in the circular 
and linear polarization, respectively (I c being the contin- 
uum intensity). 

3. SMALL DIPOLES COUNTING AND STATISTICS 

We look for small magnetic loops in the data set: co- 
herent structures that appear as a dipole in the longi- 
tudinal magnetogram (i.e., adjacent positive and nega- 
tive patches in Stokes- V maps), and a linear polarization 
patch between them (sec Figure 1), and that remain iden- 
tifiable for several (at least 2) time frames (>1 minute). 

We looked for small magnetic loops by direct visual 
inspection. A systematic search was performed in both 
field-of-views (FOVs) at all times, for these structures 
by one of the authors (MJMG), recording the position 
and evolution at different times of every single event. 
The method was validated by selecting a small area 
(10 x 10 Mm 2 ) and independently looking for such struc- 
tures by a different observer (RMS). In this control re- 
gion, the first observer found Nm — 37 loops, the sec- 
ond one Nr — 40. Nmhr — 29 of them were found by 
both observers. From those values, the actual number of 
loops N on the area may be estimated by the Laplace ra- 
tio (or Lincoln-Petersen index) fo r population es timates: 
N = NmNr/Nmhr = 51 (e.g.. lCochraj]li97l : a bet- 
ter, less biased, estimat e being N = (N m + 1)(Nr + 
l)/(N M r\R + 1) - 1 = 50 (|Chapmanl [T951'). The variance 
of this estimate is Var(TV) « (N M + 1){N R + 1)(N M - 



NmhrMNr - N M nR)/\(N M nR + l) 2 (N MnR + 2)] = 4.9 
(|Chapmanlll95l . 

These formulae have been used in the literat ure to 
estimate the size of ecological populations (e.g., iSeberl 
2002; Krcbs 2008), and the number of errors in a mes- 
sage ()Barrowl H998) .The main difficulty for applying the 
method to our case is that the objects to be counted 
(the loops as defined in the first paragraph of this sec- 
tion) may not be unambiguous. To guarantee that this 
condition is fulfilled (i.e., that both observers were iden- 
tifying and counting the same objects), both observers 
went again over all the structures found in the area by 
both, and they had to agree over all the events to be 
counted as loops. 

From this analysis we conclude that the total number 
of magnetic structures that we find in the whole dataset 
could be understimated by ^35%. We note that the 
most clear cases (those with the strongest polarization 
signals, lasting longer, and clearly isolated from neigh- 
bouring magnetic patches) were often found by both ob- 
servers in the control area. The 35% discrepancy is due 
to those events found by one observer but not the other; 
these cases correspond to the most subtle (often border- 
ing the detection limit) events. 

We identify 497 small magnetic loops emerging in 
the observed regions. Taking into account the to- 
tal time of observation and the spatial area cov- 
ered, this amounts to an emergence rate of small 
dipoles of 0.25 loops h -1 arcsec -2 , rising up to 
0.34 loops h" 1 arcsec -2 when we apply the correc- 
tion for undetected (although present in the data) 
loops given above. These values are one order of 
m agnitude larger than the previous estimation found 
bv IMartinez Gonzalez fc Bellot Rubiol (|2009f h but com- 
patible with Martinez Gonzalez (2007) and Martinez 
Gonzalez et al. (2010). There are several reasons for 
the discrepancy. First, the former study covered a rela- 
tively small area and, if the emergence is not strictly uni- 
form (as we shall discuss below), the rate can be greatly 
understimated. This fact was already p ointed out in 
IMartinez Gonzalez fc Bellot Rubiol ()2009f ). where it was 
noted that there was evidence for preferential emergence 
regions. The reason why the studies in the near-IR and 
this letter are in agreement is because, either having 
a better Zeeman sensitivity (Martinez Gonzalez et al. 
(2007) and Martinez Gonzalez et al. (2010) used the 
more sensitive- most importantly to linear polarization- 
spectral line of Fe I at 1.5 /zm) or having a better spatial 
resolution makes both studies more sensitive, inducing 
the identification of more (weaker and/or smaller) struc- 
tures in the field. 

Both opposite polarity feet and the linearly polarized 
bridge between them were tracked during the full loop 
phase. In most cases (60% of the events), the loop 
collides or merges to some degree with a neighbour- 
ing structure and we could not trace further its indi- 
vidual evolution. This is a common case because at 
our level of sensitivity and at such spatial resolution, 
most of the pixels in the FOV show circular polarization, 
and many, both linear and circular. For the remaining 
40% of the cases, the loop individual history could be 
traced beyond (and before) their complete loop phase 
(i.e., when both three polarization patches are seen si- 
multaneously). In nearly half of these cases (56%), both 
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Fig. 2. — Panel a) displays the distribution of the initial distance between footpoints (black line) and its maximum extent (blue line). 
The vertical dashed line represents the mean diameter of the footpoints at the initial time of detection (when the two footpoints arc visible 
for the first time). Panel b) displays the time evolution of the footpoints distance with respect to the initial footpoint separation. The thick 
blue line represents the average and the dashed line is the linear fit. From the slope of this fit we infer the velocity of footpoint separation 
of the small dipoles. Panel c) and d) represent the histograms of the magnetic flux and the magnetic flux density of the footpoints of the 
loops, respectively. The black lines repres ent the positiv e value s while the blue lines are the negative ones. The dashed line corresponds to 
a power law using the exponent found by Parncll ct al. (2009). Note that the proportionality coefficient is different since the fluxes have 
been obtained with two different instruments and have probably an offset. 



footpoints and linear polarization appear and/or disap- 
pear simultaneously because the structure falls below 
the detection level of the instrument (or it submerges). 
We named this population of loops as " low- lying" (see 
IMartmez Gonzalez k, Bellot Rubiq 12009) . This popula- 
tion of loops is very particular of the quiet Sun. In 42% 
of the cases, linear polarization precedes the detection 
of both footpoints and then disapp e ars be fore them too. 
IMartmez Gonzalez fc Bellot Rubiol (|2009T ) computed the 
line-of-sight velocity of the loops using the Stokes V zero- 
crossing shift and verified that all the loops having this 
very same time evolution were rising J7-loops. We have 
no reasons for thinking that the loops found in IMaX are 
a different population. However, we cannot compute the 
Stokes V velocity in the IMaX data sets analyzed in this 
paper, hence, the identification of these loops as rising 
f2-loops can be put in doubt. In just 3 instances (2%), 
we found that first the opposite polarities appear, then 
the linear polarization between them, and everthing dis- 
appears, wich could be interpreted as the emergence of 
a U-loop, or, more probably (considering the local evo- 
lution of the flow), a submergence of an f2-loop. 



Figure [T] display four examples of small loops found in 
the IMaX data. Note that, for the sake of clarity, we 
have only drawn the contours of interest, avoiding circu- 
lar and linear polarization patches adjacent to the loop 
structure. The time runs from top to bottom, the time 
cadence being irregular. The first and third columns rep- 
resent typical loops in which the linear polarization dis- 
appears at some point in the evolution of the loop while 
the footpoints stay in the photosphere (probably rising 
i7-loops). These two loops have, however, some pecu- 
liarities that differentiate them. The loop in the first 
column appears at the border of an expanding granule. 
As the granule expands, the entire loop is dragged by the 
plasma flow. The loop in the third column contains a lin- 
ear polarization signal with a gap in between. These two 
peculiarities are a consequence of the spatial resolution 
of the IMaX data that allows us to trace the dynamics 
of the linear polarization. 

The loop in the second column of Fig. [1] has an- 
other linear polarization patch with a curious dynamics. 
It seems that, somehow, the negative footpoint is dis- 
conected from the positive one (i.e., the loop breaks) and 
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Fig. 3. — Magnetograms of the two areas analyzed in this paper integrated for the whole observational period. The grey scale has been 
saturated to show just the most strongly magnetized regions (otherwise, virtually all pixels in the FOV show a det ectable signal above noise 
level) . Yellow patches show the average linear polarization signals after correcting from the bias according to ( Martinez Gonzalez ct al. 
12012 ). and therefore, they are patches with a reliable detection of the transversal component of the magnetic field. Red circles show 
the average position (between both footpoints) of the 219 (panel a) and 216 (panel b) small scale loops detected in both datasets. The 
trajectories they follow while they evolve are comparable to the size of the circles. The two large areas covered by the dotted circles in panel 
b) mark the two dead calm regions with a statistically significant lack of (detectable) magnetic loops. Note that they show low magnetic 
activity too. Arrows indicate a few hotspots where several loops appear and disappear at approximately the same position, succesively 
along an extended period of time. 



that it connects somewhere else to the right of the pos- 
itive footpoint. Of course, it is just a visual impression. 
The example in the fourth column is an example of two 
loops appearing very close in time and in space. They 
also disappear more or less at the same time, hence, one 
could think this is an evidence of a sea-serpent magnetic 
field line. Note also that the rightmost loop rotates. 

All the observed dipoles are smaller than ~ 1 Mm 
(center-to-center distance between the two opposite po- 
larity patches), becoming increasingly more abundant at 
smaller scales, with most of the observed dipolar struc- 
tures being ~ 0.4 Mm (Figure 2a). This is barely three 
times the spatial resolution limit of IMaX, which sug- 
gests that the detection and our statistics might be lim- 
ited by the instrument. The tilt angle of these dipoles 
is uniformly distributed, meaning that it does not fol- 
low the Hale's polarity law. This last resul t is co nsis- 
tent with lMartinez Gonzalez fc Bellot Rubiol (|2009t ) and 
even with the behaviour of the smallest ephemeral re- 
gions. Figure 2b displays the distance between foot- 
points with respect to the separation between them at 
the initial time. Therefore, positive values mean that 
the footpoints separate with time and negative ones in- 
dicate that the footpoints approach each other. In aver- 
age, the distance grows linearly with time with a velocity 
of Vd — 0.9 km s _1 , comparable to typical granular val- 
ues. This indicates that the loops passively follow the 
granular flows, as expecte d from weak magnetic features 
(jManso Sainz et al.ll20Toh . 

The magnetic properties of these small dipoles are rep- 
resented in Fig. 2c and 2d. They are obtaining invert- 



ing the data in the weak field approximation following 
lMartmez~G onzalc z et al.l (|2012D . The magnetic flux has 
been computed in the area containing the observed sig- 
nal. The frontier has been defined by eye (as an iso- 
contour of magnetic flux density) and hence is slightly 
different for the different structures. The magnetic flux 
density is the mean value of the magnetic flux densities in 
this same region. The magnetic flux of the loops can be 
explained wit h a po wer law using the exponent found by 
iParnell et ail (|2009l ). This means that the population of 
small dipoles follows the population of magnetic fields in 
the quiet Sun. But looking at the histogram of magnetic 
flux densities, the loops are located at the end tail of the 
histogram; it is mostly in the range 10-20 Mx cm -1 — i.e., 
10-20 G if the magnetic field were uniformly distributed 
a nd volume filling. This is comp atible with the results 
of lMartinez Gonzalez et ail (|2010D who state that, in the 
quiet Sun, the larger the signal the larger the degree of 
organization of magnetic fields. 

4. SPATIAL DISTRIBUTION 

Although small scale loops are found all over the ob- 
served areas, their spatial distribution does not seem to 
be completely uniform (Fig. 3). It can be observed that, 
at some locations, loops appear repeatedly and succe- 
sively as in bursts, forming clusters, a behaviour that has 
been noticed before by Martinez Gonzalez & Bellot Ru- 
bio (2009), who pointed out that, often, the appearence 
of loops made it more likely that new ones were later de- 
tected nearby. On the other hand, extended areas seem 
to be noticeably empty of such events, as if voids ap- 



Dead calm areas in the quiet Sun 



5 




A (Mm') 

Fig. 4. — Number of voids per unit area in the two datasets. 
Dark blue and cyan histograms correspond to the data shown in 
Figures 3a and 3b, respectively. Solid line shows N(A) for a Poisson 
distribution of points with n = 0.22 Mm -2 (we use a single value 
since we observe almost the same number of loops in both data 
sets, i. e., 248 and 249 loops). The inset window represents the 
probability of finding an area > A. The x-axis is the area A in 
Mm 2 . 

peared in the distribution. However, this may be deceiv- 
ing since voids are also formed even in strictly un iform 
distributions of points (|Betancort-Riidll990l fl99l . 

A quantitative analysis to determine if these voids are 
statistically significant was performed. This requires on 
the first place, an unambiguous definition of 'void', a 
non-t rivial task in itself (see e.g . . iKauffmann fc Fairalll 
fl99lt iTikhonov fc Karachentsevl 120061 and discussions 
therein). We adopted the simplest definition here and 
considered only voids of circular shape: the largest empty 
circle that can be fitted in a given region of the point field 
— equivalently, an empty circle limited by three points of 
the distribution^]. 

If the loops appeared uniformly on the solar surface, 
then the probability of finding a void of area between A 
and A + dA within the field of view (square surface of 
area L 2 ) would be (see Appendix): 

P(A)dA = —3 — {L r y^e~ nA ndA, (1) 

-Avoids A 

where A = 7rr 2 , n is the surface density of points, and 
Avoids i s the total number of voids espected in the FOV, 
which is given by Equation (Al). For the two data sets 
studied here L = 31 Mm (which is slightly smaller than 
the nominal FOV because we have excluded the apodized 
exterior area), and n — 0.22 Mm~ 2 . Note that we only 
need a single value of n since the number of loops de- 
tected in both data sets are very similar (i.e., 248 and 
249 events). 

Figure 4 shows the number of voids per unit area in 
both datasets and for the corresponding Poisson distri- 

1 This is certainly an overabundant definition: several overlap- 
ping circles may be found covering what we intuitively consider 
as a single "void". Algorithms might be devised to merge circles 
and to find a definition closer the intutitive meaning ( Gaite] 120051 ; 
IColber g ct al. 2008;). This "overcounting" is, however, of no im- 
portance for our calculating the probability of finding a void larger 
than a given size (roughly, all voids are overcounted equally), and 
we will use this much more simple approach which avoids numerical 
technicalities. 



bution. Except for the smallest areas, the distribution of 
large circular voids in the first dataset is not significantly 
different from the uniform one. On the second data set, 
however, apart from the overabundance of small scale 
voids, large circular voids seem to be significative ly more 
abundant than a strictly uniform distribution would sug- 
gest. Actually, for the parameters found for these obser- 
vations, the probability of finding a circle with an area 
lager than A = 35 Mm 2 (equivalently, diameter larger 
than 3.3 Mm) is very low (3 x 10 -4 ; see the inset plot 
in Fig. 4). We conclude, therefore, that there is statis- 
tical evidence for the two voids marked in Figure 3a to 
be real and not due to chance. Moreover, the oberabun- 
dance of small voids is interpreted in terms of a clumpy 
structuring (see how loops appear in clumps, like a gurgle 
phenomena, in Fig 3). 

In order to relate the distribution of loops (and the 
voids) to the global magnetism in the observed area, Fig. 
3 shows the integrated longitudinal (in black and white) 
and transverse (in yellow) magnetograms for the two ob- 
servations. When dealing with the linear polarization, 
one has to remind that it is a biased estimator of the 
transversal field component. In the plot, this bias is sta- 
tistically partially removed as follows: we compute the 
bias f or a percentile 95 when the ob servations are pure 
noise ([Martinez Gonzalez et al.ll2012L see) (note that this 
bias value depend on each pixel, i.e., on the actual in- 
tensity profile). This means that the "real" transverse 
magnetic will be below this bias value with a probability 
of 95%. We have decided to put all the values smaller 
than this bias to 0. Figure 3 shows only the statistically 
significant patches of linear polarization appearing at all 
the observed times. The positions of the loops do not 
seem to be clearly correlated with the longitudinal mag- 
netogram, but the voids encircled by the loops show less 
magnetic activity than other areas in the FOV. Consid- 
ering the correlation with linear polarization, it suggests 
that most of the linear polarization signals that are de- 
tected are associated to loop structures embeded in the 
formation region. 

5. DISCUSSION 

It is known that even in very quiet areas of the Sun 
magnetic fields may organize naturally forming loops 
at granular scales. In this study we extended this ob- 
servation to the smallest spatial scales observable (100- 
1000 km), finding an increasing number of loops at 
smaller scales up to the resolution limit. This finding 
suggests that the organization of magnetic fields might 
continue beyond that limit. We cannot reconstruct the 
complete magnetic field topology because 1) the finite 
spatial resolution of our observations is (perhaps inher- 
ently) above the organization scale of the magnetic fields, 
and 2) we lack linear polarimetric sensitivity, which gives 
us only fragmentary information on the transversal (hor- 
izontal) component of the magnetic fields. Due to these 
limitations the loop structures that we observe are bi- 
ased towards relatively large and relatively strong with 
respect to the magnetic flux density in the neighbouring 
areas. We found evidence that the loops thus detected 
are not randomly distributed on the solar surface, but 
rather that they may appear in bursts, and that they are 
noticiable absent from extended areas which are, also, 
only weakly magnetized. 
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It is not yet clear what is the nature of the magnetic 
fields in the quiet Sun — what are the fundamental phys- 
ical mechanisms involved in their generation and evolu- 
tion. The presence of these dead calm areas in the quiet 
Sun (and small scale loops hotspots) represent an impor- 
tant constrain on the origin of magnetic fields in the very 
quiet Sun and on the dominant dynamic and magnetic 
mechanisms taking place there. 

It is thought that the magnetism of the quiet Sun can 
be the result of the emergence of und erlying organized 
magnetic fields (|Moreno-Insert isJ l2012fl or the dra gging 
of the overlying canopy fields (jPietarila et al.l 120111 ). It 
would then be necessary to understand why there are 
emergence hotspots and dead calm areas. Another pos- 
sibility is that they are just recycling of the decay of 
active regions as they diffuse and migrate to the poles. 
But it seems unlikely that such random walking would 
lead to the kind of organized structures and to the spatial 
patterns reported here. It is also possible that they are 
linked to som e type of dynam o action taking place in the 
solar surface (|Cattaneolll999[) . It seems now clear that 
coherent velocity p atterns are a requisite for dynamo ac- 
tion to take place ([Tobias fc Cattaneol 120081 ). The most 
obvious coherent velocity pattern in the solar surface is 
granulation. If this velocity pattern is involved in some 
dynamo action, it is reasonable that it forms coherent 
magnetic structures (such as the loops), although these 
does not need to be organized at the same granular scales; 



it could well be that they form intermittent patterns as 
the o nes observed here. Ac tually, theoretical considera- 
tions dChertkov et al.lll999l ) and laboratory experiments 
(jRavelet et al.ll2008h support the idea that the onset of 
turbulent dynamo action may be highly intermittent and 
bursty. Finally, it could just be that these small scale 
loops represent the far tail of a continous range of struc- 
tures from a global dynamo, just lying at the other end 
from sunspots. Their spatial statistics would then reflect 
the velocity patterns on the last (shallowest) layers of 
magnetic field emergence. 

Future models that we construct to understand the 
generation of magnetic fields in the very quiet Sun have 
to explain the spatio-temporal coherences that we report. 
Further work is needed to extend these results in larger 
areas of the Sun and along the solar cycle. 
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APPENDIX 
DERIVATION OF EQUATION (1) 

We derive Equation (1) adapting to our case some arguments of the strategy bv iGaitd ((2009). The probability that 
three points i — 1, 2, 3 distributed at random (uniformly) on an square with area A t = L 2 , have coordinates between 
Xi and Xi + dxi, yi and yi + dyi (0 < Xi < L and < yi < L), is (l/A 3 )d 3 Xid 3 yi. With the change of variables 
Xi = x c + r cos 8i, yi = y c + r sin 8i, we may express the probability for the center of circumcircle of the three points to 
lie between x c and x c + dx c , y c and y c + dy c , its radius between r and r + dr, and the azimutal angles of each point 
between 8i and 8i + d9i, then, as dx c dy c drd 3 8i\ J(r, 8i)\/A 3 , where J(r,8i) = r 3 (cos 8\ (sin 82 — sin #3) + cos (^(sin 83 — 
sin 8 1 ) + cos #3 (sin #1 — sin #2)) is the determinant of the Jacobian matrix. 

The probability of one such a circle to lie between the bounds of the large square and its radius to be between r and 
r + dr is P m (r)dr = 24tt(L — r) 2 r 3 dr/A 3 , where x c and y c have been integrated betwen r and L — r, and 8i between 
and 2-7T, taking into account that the integral of the angular part of |J(r, 8i)\ is 24ir 2 . Alternatively, the probability 
density of the area A = ixr 2 of such a circle is P\ n (A)dA = 12(L — r) 2 AdAj A\. 

On the other hand, in a homogeneous Poisson field with density n = N t /A t (N t points randomly distributed over an 
area A t ) , the probability of finding k points in a region of area A is P k (A) = (nA) k /k\e- nA . Therefore, the probability 
of the circumcircle of three points to be void is 12(L — r) 2 AP a (A)dA/ A 3 . 

Finally, N t points determine N t (N t - l)(N t - 2)/3! « N 3 /3\ (N t > 1) different triplets (hence, possible circles), and 
the number of void circles with areas between A and A + dA within the square bounds is 2(L - r) 2 Ae~ nA n 3 dA. 

The total number of void circles is: 

r Nm 2 

N voids = 2(L-r) 2 Ne- N ndN=-[2{l + 2N m )-e- N ™(2 + Nl)-3y/N^erf[y/N^}}, (Al) 

Jo ^ 

where N rn = mrL 2 /4 and erf(x) is the error function (Ab ramowitz fc Stegunl[T972[ ). The probability of voids of area 
A with our constrains is then the number of voids of area A divided by the total number: 

P(A)dA = —?—(L - r ) 2 ^^e- nA ndA. (A2) 

-Avoids A 

For a large total area A t (L,A t ,N m — > 00 at constant n), then P(A)dA — ^ nA ^ e~ nA dA, which coincides with the 
analysis of iPolitzer fc Preskilll (|1986f ). 
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